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ABSTRACT 

Parameter plane techniques of analysis and design are applied to 
variable parameter linear electric networks. Parameter plane theory is 
extended in several areas. 

A theorem permits the constant zeta and omega loci to be obtained in 
non-parametric form. 

The class of network functions to which parameter plane techniques 
may be applied is broadened by permitting parameters to appear non- 
linearly in the polynomial coefficients. 

Algebraic design methods are presented which permit the solution for 
parameter values which will simultaneously meet many different types of 
specifications. 

Pole-zero sensitivity formulas are derived from the non-parametric 
equations of the constant zeta and omega loci. 

Equations for curves of constant bandwidth and Q are derived which 
permit these characteristics to be displayed on the parameter plane. 

Examples are presented which demonstrate the application of these 


techniques to several common networks. 
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bi Introduction 

Linear active and passive electric networks frequently have one or 
more variable parameters. In certain instances network elements are in- 
tentionally varied so as to modify the characteristics of the network. 
In other instances the network has no intentionally variable parameters, 
but due to component tolerances, aging, heating etc., variations do oc- 
cur. If it is assumed that these variations occur sufficiently slowly 
so that the network may be considered time invariant, and if it is de- 
sired to predict or control the effects of such variations on the charac- 
teristics of the network, then the methods to be described can be applied. 

Parameter plane techniques were first introduced as a technique for 
analysis and synthesis of feedback control systems. They were used to 
investigate the effects upon the roots of the control system character- 
istic equation as two parameters within the system were varied. 

In the early 1950's Mitrovic (1) considered system parameter varia- 


tions which directly affected two of the coefficients of the system 


characteristic equation. For example: 


ee te iesete te amar +a a + B.s + B. = 0 (Tha) 


n 
ae 3 2 1 0 


n-l 


where: as are constant and real 


Bo and By are variable and real 


Mitrovic's original work resulted in curves of constant zeta, omega, and 


sigma which he drew upon the B “By plane and which specified the roots of 


0 


Equation (1-1) for any choice of By and B,- 


In 1964, Siljak (2) extended Mitrovic's work to include cases where 


variable system parameters, say o and 6 appeared linearly in one or 





more of the coefficients of the characteristic equation. Specifically 


Siljak investigated the case where: 


~s ais = 0 (1-2) 
L=O 
where: 
MP ED eta etd. 
i i i 
and: 


b,,; @., d. are constant and real 


«x, @® are variable and real 
Siljak, like Mitrovic, obtained curves which were plotted upon the X-8 
plane and which specified the roots of Equation (1-2) for any choice of % 
and 6 . It was Siljak who first referred to such graphs as "parameter 
planes". 

The parameter plane provides information about the roots of the 
characteristic equation similar to that provided by root loci, but unlike 
root locus techniques, which investigate the effects of varying one system 
parameter, parameter plane techniques are inherently best suited for in- 


vestigating the effects of varying two system parameters. 





Zz. Basic Parameter Plane Techniques 

Parameter Plane techniques evolved from a study of the characteristic 
equation of control systems and of how the roots of such equations were 
affected as two parameters of the control system were simultaneously varied. 
More basically however, the parameter plane is a method for investigating 
the roots of variable coefficient polynomials, and hence is applicable in 
a much wider sphere than that of control systems alone. This chapter will 
review the parameter plane techniques as presented by Siljak (2) and will 
point out some of the reasons that have limited the scope of applicability 
of these techniques. Subsequent chapters will present new work which 
overcomes some of these limitations. 


Let us consider a polynomial equation: 


an k 
F(s) = Y as = 0 (2-1) 
kK =0 


where the coefficients a, (k = 0, 1, ...,n) are real and s is the 
complex frequency variable s = @ + je = - Fw+jo4 1 -3? 
If one forms powers of s it can be shown (2) that the following relation 


holds: 
sew n(-7 +5 F1-3? UC-7)) (2-2) 


where: 


k 
T, ¢- % ) = (-1) T, ( 3) 
(2=2) 


U.(-F = CD CZ) 


The UL and Ty are Chebyshev functions of the first and second kind 


respectively and are given by the following recursion relations (2): 


“isso 
ge ee = et — 
—— . oz ———— as 60 

ci Om — + = eee 61 

> -— =a . ————— 

ee — ——_— <r —_ a 5 

— ; ne el 

— —) <-> € => _ 
=~ —— =< & 


C7 ® s» _ /—~ 


- © —_ ~~ 








Be (Op at) FT OZ ES, 0 
(2-4) 


U4 CF) 2 UCT +Y CF) = 0 
where: 


To Fd = bs TF 27 UF) =0,0 0%) = 1 


Tables of these Chebyshev functions are included in Appendix III. 
Substituting Equation (2-2) into Equation (2-1) and setting the real 


and imaginary parts of the resulting equation equal to zero independently we 


5 a, WT, ( -4) 


K=O 


obtain: 


il 
© 


(2-5) 


il 
=) 


P 
k 
a, WU, ( -3) (2-6) 
K=0 
Employing Equations (2-3) and (2-4) we obtain these expressions in 


terms of the tabulated Chebyshev functions: 


nN 
) (-1) “ay oD, CF ) = 0 (2-7) 
K=0 

kk 
3 (-1)"a,@°U.( 5 ) = 0 (2-8) 
K=0 


Let us consider, as Siljak did, that the coefficients a, of Equation 


k 


(2-1) are linear functions of the variable parameters & and @ as follows: 


a. = bo + c. B + qd. (2-9) 


where: i 


b dy. (k = 0, 1, ...,n) are constant and real 


» Cy, 
k° ok 


Blue € are variable and real 


es Lo, r 
1 a (ee 





Substituting Equation (2-9) into Equations (2-7) and (2-8) we obtain: 


BX +C,@ =-D, (2-10) 
BX wr =~ % 
where: 
B, = )) (-1), ou, Bare S (-1) bw “U, 
=O K=0 
k , kk 
oF (-1)%c, WU, Con S (-1)"e, WU, (2-11) 
K=O K:0 
k \ k, .k 
D, = a (-1) “a, wD, D, = S (-1) “a, ©", 
K=0 K=O 


Equations (2-10) are two linear equations in the variables %™ and 8 F 


Solving these equations we obtain: 


of = © P9°SoPy nae as i tae (2-12) 
B,C,-B,C, B,C, -B,C, 


Equations (2-12) give o& and éG in terms of 3 and 4 . Thus if we 
fix A on and vary W over the range 0 < W< © we will ob- 
tain a curve in the &- 6 plane which specifies the relation between & 
and @ necessary to cause a pair of complex roots of Equation (2-1) to 
have a damping ratio "4 equal to Zo: Curves resulting from fixed 
values of are referred to as constant zeta loci. If we fix W= Wo 
and vary 4 over the range -l < Z < +1 we will obtain another 
curve in the o&- 8 plane which specifies the relation between & and /3 
necessary to cause a pair of complex roots of Equation (2-1) to have an 


undamped natural frequency W= JW) Curves drawn for fixed values of 


0° 


W are referred to as constant omega curves. 





If we substitute s = @ in Equation (2-1), where sigma refers to 


values of s along the real axis in the s-plane, then we obtain: 


nN 
»y atom = © (2313) 
K=Q 
If we now substitute for ay the linear relation (2-9) we obtain: 
n n 
>) bv + BY yo" + Sy aot = 0 (2-14) 
K=O =O K=O 


For fixed values of sigma,Equation (2-13) represents a straight line in 
the & - B plane. This straight line specifies the relation between A 
and B necessary for Equation (2-1) to have a real root at s = g@ . Such 
lines drawn upon the &-f plane are referred to as constant sigma loci. 
If the reader is interested in the mapping which occurs between the 
s-plane and the parameter plane or in the properties of this mapping then 


the author refers him to reference (2) where this is discussed in detail. 


Example I: 


Let us consider that we have constructed a parameter plane for some 
third order polynomial and that this appears as shown in Figure l. 

Let us consider that parameter values & and 8 are chosen to be & 0 
and B 0 respectively, and that at these coordinates on the parameter 


plane it happens that the 0&7 x4 and 4), loci intersect. We may 


ihe 
thus assert that for this choice of parameter values the polynomial has 


the roots: 


s = v, 


ae Bi Ds Be z,7)7 


Let us now presume that we wish to change the damping ratio to oa: 








es) 
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Fig. 1. Example of a Parameter for some 3rd Order Polynomial 


but wish to leave the undamped natural frequency unchanged at W 7° 


We therefore find the intersection of the z 0 and Ud. loci, Ine thas 


example this occurs at KH = & 1? 3 = B Adjusting the para- 


1 


meters © and B to these values will guarantee complex roots at s = 


| , 2. % aap 
: Zo @ , + j Ww, (- er The real root which was ats= @ 4 


has moved to a new location between wT, end ©” 3° 


pine 
Thus we see that the parameterfis a general tool of both analysis 


and synthesis in two parameter systems. 





O Scope of this Work 

This section outlines the general content of this manuscript and 
introduces new work which overcomes some of the limitations of parameter 
plane theory which have heretofore prevented its application to electric 


networks. 


General Content 


In this work the author investigates the application of parameter 





plane techniques to variable parameter electric networks. The objective 
of this work is os provide means of analysis which will show how the net- 
work characteristics change as the parameters are varied and to provide 
means of design which will permit values of the adjustable parameters to 
be chosen so as to cause the network to have desired characteristics. To | 
this end both graphical and analytic techniques are presented. The net- 
work characteristics which can either be determined or controlled are 
pole and zero locations, frequency response, Q and root sensitivity. 

Section 4 through 8 develope the mathematical theory which related 
the network characteristics to the variable parameters. The mathematical 


theory presented in these sections forms the basis for both the graphical 


‘ ‘ 
e 
@ 


and analytic techniques. 

Sections 9, 10, and 11 contain examples of graphical analysis and de- 
Sign of selected networks which have historically been of considerable 
interest. These examples provide not only an illustration of the poten- 
tial of the parameter plane techniques, but also provide normalized results 
which will aid designers in employing these networks. 

Section 12 considers several other networks of general interest. 


These networks are considered in enough detail to show the applicability 





of parameter plane techniques to their analysis and design. 

Section 13 presents an algebraic method of design which employs the 
parameter plane equations developed in earlier sections. These algebraic 
methods of design can be applied to any problem which may be graphically 
designed and additionally to the much larger class of problem which have 
more cae two parameters. There is theoretically no limitation on the 
number of parameters which may be simultaneously considered when using 
these algebraic techniques of design. 

Section 14 concludes the manuscript with the author's comments on 
the work and with his recommendations for additional work in the area of 


network analysis and design by parameter plane techniques. 


Limitations of Earlier Parameter Plane Theory and Introduction to New Work 


Parameter plane theory as it has heretofore existed has several 
limitations which have restricted its application to electric network 
analysis and design. The remainder of this section discusses these limita- 


tions and introduces new work aimed at overcoming them. 


Form of the Coefficients 
Consider the restrictions which have been placed on the polynomial 
coefficients - namely that the variable parameters must appear linearly in 


the coefficients: 


ay = b & + c 6 & qd. 


This restriction was not too severe when parameter plane techniques were 
applied primarily to feedback control systems since many useful systems 
have characteristic equation coefficients of this form. Electric networks 


on the other hand, seldom have coefficients of this type. Consider for 


9 








7 





example the Parallel-T network shown in Figure 2: 


m7 


AC QR = 


Fig. 2. Symmetrical Parallel-Tee Network with Variable Elements 


This network has the transfer function: 


E (Pp) = X68 P -++ 2 B se + 2 8 p + |] (3-2) 


Cw Oo ee Oe eee 





E, (p) XB eek (2B +24 +X )p? +(28+% +2) +1 


where; p = RCs 

Observe that Equation (3-2) contains terms (6 as well as the linear 
terms permitted by Equation (3-1). Thus we are unable at this point, due 
to the linear requirements placed on the coefficients, to determine by 
parameter plane techniques how the singularities of (3-2) are altered as 
the parameters o& and A - are varied. New work to be presented in 


Section 4 will overcome this difficulty. 


Constant Zeta and Omega Loci 


As derived by Siljak (2), the constant zeta and omega loci are 
presented in parametric form as functions of the independent variables 


omega and zeta respectively. This may be seen in Equation (2-12) where 


10 








the constant zeta loci are obtained by fixing zeta at selected values and 
generating the loci by varying the independent variable omega over the 
range 0 to infinity. Similarily, the constant omega loci are obtained 
by fixing omega at selected values and then varying zeta as an independent 
variable over the range -1l to +l. 

Clearly it is desirable to obtain the equations of these loci in non- 
parametric form. By achieving this we simplify the plotting of these loci, 
and more important we obtain these loci as families of plane algebraic 
curves in the & - ® plane. This will permit the mathematical procedures 
for the investigation of plane algebraic curves to be applied when analyz- 
ing the behavior of these loci. This is an important step if procedures 
for sketching the constant zeta and omega loci are to be developed. In 
Section 5, procedures are presented for reducing these equations to non- 


parametric form. 


Frequency Response Characteristics 


Since frequency response characteristics are often used in specifying 
the behavior of an electric network it is essential that any method which 
purports to be a general tool for analysis and design of such a network 
include procedures for determining and controlling the essential features 
of its frequency response. Choe (3) developed formulas for displaying 
curves of constant bandwidth upon the Mitrovic plane, but these techniques 
were not applicable to the parameter plane. In Section 6 formulas are 
developed for displaying curves of constant bandwidth upon the parameter 
plane. These formulas are not only applicable to electric networks, but 
to any rational system function which can be displayed upon the parameter 


plane. 








Constant Q Loci 


Second order network transfer functions can be characterized by their 
Q, much in the same way that the characteristic equation of a second order 
servomechanism is characterized by its damping ratio. Procedures are 
developed in Section 7 for displaying loci of constant Q upon the parameter 
plane. These specify the value of Q which results for any choice of network 
parameters or alternately specify in what fashion the network parameters 


must be changed in order to obtain some desired value of Q. 


Root sensitivity to Parameter Variations 


The sensitivity of the roots of a polynomial to changes in parameters 
contained in its coefficients can be catagorized as macroscopic or microscopic. 
Macroscopic sensitivity refers to the sensitivity of the roots to large changes 
in the parameters. Microscopic sensitivity refers to infinitesmal changes in 
the parameters. The constant zeta and omega loci drawn upon the parameter 
plane permit macroscopic sensitivity to be determined by inspection. Micro- 
scopic sensitivity may be determined from formulas derived by Kokotovic and 
Siljak (5) from the parameter plane equations or may be computed directly 
from the non-parametric equations of the constant zeta and omega loci which 
are developed in Section 5. Section 8 describes how both macroscopic and 


microscopic root sensitivity to parameter variations may be determined. 


Examples of Parameter Plane Techniques Applied to Networks 


Sections 9 - 11 contain detailed examples of parameter plane techniques 
applied to specific networks. These examples were carefully chosen in 
order to both fully demonstrate the current state of parameter plane tech- 
niques, and at the same time to produce results useful to persons desiring 


to employ these networks. The results which are presented are general and 


12 











are presented in normalized form. 
Section 12 contains examples of several networks with variable para- 
meters which were studied only so far as to show the applicability of 


U 


parameter plane techniques to their analysis and design. 


tS 








4, Parameter Plane Equations for Coefficients of the Form: 


a, = bX +e, 8 +d «6 + £) 
In this section equations will be obtained which will permit para- 


meter planes to be constructed for polynomials whose coefficients are of 


the form: 


= SPAS Tc ae ie oe xB +f (4-1) 


ak k 


where; b d. are constant and real 


b> “K? WK 


Oo, 8 are variable and real 


Consider the polynomial equation: 


Nn 
me) i as. at (4-2) 
K=0 
Let the coefficients aL be of the form specified by Equation (4-1). 
Let: 

k= k 2 

s = WIT,(-F +i V1- 7% U.(- 3d} 4-3) 
where: T, and U, are Chebyshev functions discussed in Section 2. 


k k 
Substituting (4-3) into (4-32) we obtain: 


n n 
F(s) =) aw -F Vt GfI-F? Y awh yg ) 
K=0 K=O 


(4-4) 
Equating the real and imaginary parts of Equation (4-4) to zero independent - 


ly we obtain: 


Y\ 
> 4 WS T(- 5) = 0 (4-5) 
K=0 
S aw U(-% ) = 0 (4-6) 
K=0 


14 








but from Equations (2-3) and (2-4): 
T = 1 bs 


u-3)= CD ucy) (4-7) 
k 
Teal % ) 6GD" YD 
Substituting Equations (4-7) into Equations (4-5) and (4-6) we obtain: 
< k k 
» Cire. 12 ) 
K=O 
(-1)"a Ww uC e ) = 0 (4-9) 


K=90 
Substituting for the coefficients a 


il 
© 


(4-8) 


k in Equations (4-8) and (4-9) the 


form specified by (4-1) we obtain: 


n n 
KS CHO WNC e+ 8 > (-1)“c, wt, ¥ ) + 
K=0 K=0 


XB sj (-1) “ay wt (3) + S (-1) "€, wt (F) = 0 (4-10) 
K=0 K=0 
k 
XS) G1)", wu ) + ey (-1)“e, wit, ) + 
K=0 
4a» (-1) “a, wu, ( rf ) + Ss (-1)"£, uw) U, ( z ) = 0 (4-11) 
K=0 
which may be written: 
BUX +¢, +D, X86 + F) = Q (4-12) 
B,X +, B + D,, (8B + ee 0 (4-13) 
where: 
‘ k k a k k 
B, = (-1)"b, @ “TF ) B, = x (-1)"b, w U,( F ) 
K=O 


(-1)"c, wt, (% ) Go. = bi (-1)e, a U,( 5 ) 
K=O (4-14) 


OQ 
A 

it 
Ms 2 


x 
tt 
o 
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(-1)"a, UF ) 


és 


k k 
D, = : (-1)"d) w T C4 ) D, 
2, k k k k 
F, » (-1)"£, w “T,¢ % ) F, = (-1)"£, WUC F ) 
K30 .@) 
Eliminating o from Equations (4-12) and (4-13) we obtain: 


x 
és 


fare =, ft Oem Ee (4-15) 
BL + DLA B., +D, B 


Which may be written: 
a ; 
Acp ®° + ‘Agp + Acs) 6 ‘Agu? “al 
F 3 


Where: 


3 ee - 1 ee 
Dep > \rp = ! — 


Equation (4-16) has the solutions: 


2 
Bs “(Ag + Ac) + (App + Aop) > “Lop Sep (4-17) 


“AQ 
Let us define: 


2 2 
“O° CIN ge 6 IN Se a 


This permits us to write Equation (4-17) as: 


ee 
hale: [Xe [ -(Agy + Acg) + 8 } Sia 


Proceeding in a similar fashion we may eliminate 6 from Equations (4-12) 


and (4-13) to obtain: 


Fee ms = 8 (4-20) 
C, + Dix C, + Dix 


As before, we may rewrite this as 
2 
Arp % + (Age t App) & + Age = 9 Geen 
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Solving Equation (4-21) for we obtain: 


1 
a = Da [ CAS FUN 2.) oh. J is 
Where we have defined: 
2 2 
Ry * (Age +Agp)” - ‘A gp Orc ei 
ill h h Re = : 
We will now show that 2 = Rg 


From Equation (4-23) we write: 


2 
Ry = (B,C,-BAC, + Dy 


Expanding Equation (4-24) we obtain: 


ae 2 2 Deseo 2.2 22 
Ree B,C, +eB, C) + DOF, + DyF, + 4 Bcaa Eee 4 BAC.D,F, seg B,B,C,C,, 
2 B,C,D,F, - 2 BiGsD Fs = 2 BC, DF, = 2 B,C,D,F. - 
) D, DoF iF. (4- 


Expanding Equation (4-18) for Re we obtain the same expression as in 


the right hand side of Equation (4-25). Thus: 


Thus we may state from Equations (4-19) and (4-22); 





1 
XK = ae (ae +f.) + R] ( 
6 = > — I-A, +A) +8) (4- 
aT FD CB 
Noting that A\..A =-Aog we may rewrite Equations (4-27) and (4-28) 
as: | 
A= oe [-(AW+ A y+ cn = [ i 
2 a. FD BC , 1 
| 0 
1 a m 
- [ -( #N/\ ..) + (-1) 8] jas! (4- 
p Won Arp BC 


17 


Z 
-B,C, + D Fy - DF.) - 4(F,C. - F,C,) (B,D, - D,8.) (4- 


R | (4- 


22) 


23) 


24) 


25) 


26) 


27) 


28) 


29) 


30) 





Equations (4-29) and (4-30) define four possible solution fit. depending 
on the values of k and m. 

chia now determine which of these possible solutions are actually 
solutions to Equations (4-12) and (4-13). Multiplying Equation (4-12) by 


D, and Equation (4-13) by -D, and adding we obtain: 


2 1 


enc *aee * Aes ° CooL) 


Equations (4-31) must be satisfied by any of the pairs (4-30) which are to 
be solutions to Equations (4-12) and (4-13). Substituting Equations (4-29) 


and (4-30) into Equation (4-31) we obtain: 


k m a 
: Ax : Axe A CG) R - Am + Asc ft CY R +A, =0 





(4-32) 
Which after simplification becomes: 
k m 
(-1)° + (-l) =0 (4-33) 
In order for (4-33) to be satisfied it is necessary that 
k - m (4-34) 
Thus we may rewrite Equations (4-29) and (4-30) as: 
X= se OL -(Aep t Age) + C'S! (4-35) 
LA sp 
1 k+1 
So sa  - (Ae Ae ee Gel (4-36) 
k 2/\cp FD BC 


Where k = 0 and k = 1 generate the two solution pairs of Equations (4-12) 
and (4-13). 

Equations (4-35) and (4-36) may be used to program a digital computer 
to generate the constant zeta and omega loci. The basic steps to be pro- 


grammed are: 
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Constant Zeta Loci 

1. Fix ] = ae This makes Equations (4-14) functions only of 
the independent variable omega. 

ae Increment omega and for each value of omega compute the solu- 
tion pairs ( 1? B,)3 k = 0, 1 using Equations (4-35) and (4-36). 

oa Connect all solutions (& 0? B 9) in the order of increasing 
omega. This generates one branch of the constant zeta locus for z = Pe: 

4. Connect all solutions (# tr a DP in the order of increasing 
omega. This generates the other branch of the constant zeta locus for 

$= 7. 
ye Repeat steps 1-4 for Z = Sts ie ”) i+2? etc. until as 


many constant zeta loci are obtained as desired. 


Constant Omega Loci 

l. Fix M= 4g) « Thus makes Equations (4-14) functions only of 
the independent variable zeta. 

a Increment zeta and for each value of zeta compute the solution 
pairs ( & le? & 3 k = 0, 1 using Equations (4-35) and (4-36). 

3. Connect all solutions ( 0° fS' 9) in the order of increasing 
zeta. This generates one branch of the constant omega locus for W = W) .. 

4. Connect all solutions ( & & D in the order of increasing 
zeta. This generates the other branch of the constant omega locus for 

W= W as 


5. Repeat steps 1-4 above for 4) = W) «449 W= yw jt2? etc. until 


as many constant omega loci are obtained as desired. 
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Dx Non-Parametric Equations for the Constant Zeta and Omega Loci. 


In this section procedures are presented which show how the equations 
of the constant zeta and omega loci, regardless of the form of the co- 
efficients a,, may be reduced to non-parametric form. 

As developed in Section 4, the parameter plane equations for co- 


efficients a. = bX +c, 3 +d) ff + £, take the form: 


x =A 7,W) (5-1) 


8 =F.) 


The constant zeta and omega loci are obtained from these equations 


by fixing one of the variables, zeta or omega, and considering the other 
variable as independent. Varying this independent variable over the range 
of interest generates the desired locus. Our objective in this section 
is to obtain the equations of the constant zeta or omega loci from (5-1) 
by eliminating omega or zeta, respectively. Achieving this will produce 


the non-parametric equations of the desired loci: 


20% , ke (5-2) 
WORX »B,W)=0 (5-3) 
Two theorems, basic to this development, are proved in Appendices I 


and II. These are: 


Theorem I 


The Chebyshev functions U,.¢ s ) obey the following relationship: 
UCF )U,9O5) - 9 C4 Gaz) = — tee (5-4) 


for all integers i and j. 
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Theorem II 


ce: 
k k 2 
S (-1) a, W U6 F ) = 0 
and if: - 
k k 
¥ (-1) ay aU ( $ ) = 0 
Then: ai 
k k 
S GPa me Wa (5-5) 
K=O 
for all integers k and j. 
Equations of the Constant Zeta Loci 
Using Theorem II for the cases j = 0, j = n we obtain: 
f(W) = (-1)"a, WUC ) = 0 (5-6) 
p U6 
K=0 
g(W) = (-1)"a, wu, F) = 0 (5-7) 
k k-n 
K=0 


Equations (5-6) and (5-7) can be considered to be polynomial equations of 
order n in W . We wish to eliminate W from these equations in order to 
obtain Z(a,, x) = 0, which is the non-parametric representation of the 
constant zeta loci. This may be accomplished by employing Sylvester's 
Dialytic Method of Elimination (6). The basic steps are shown below for 
two polynomial equations of differing order in 4) , from which 4@) is 

to be eliminated. 


Given the polynomial equations: 


f( W)=£f W)” + £ w + .eeeten0 (5-8) 


n n=1 
g( w& ) g wu) + g W + es + nO (5-9) 
n-l 0 
The quantity y) may be eliminated from Equations (5-8) and (5-9) by the 


73k 





Re RE a A TEES epee} 


following secuence of operations; 


n-2 


ie Multiply (5-8) by Mes a) Peer tae ees yon 8 vag ee 


Zs Multiply (5-9) ee ae wom See... 1 Teutensie, 

3. Form the determinant of the coefficient matrix of the 
resulting set of (mm) equations which are linear and 
homogeneous in the (nim) quantities eu). pain aes 

25, a , Te 

4. This determinant, when forced equal to zero, establishes the 

relation among the f.. and g, necessary if Equations (5-8) 


and (5-9) are to have common roots. 


The resulting determinant has the form: 


= C 24 mi * e ° ° ° e ° O 

“nm ° ° ° is 7 aC 0 

O ie ° Py e 25 a Ty O e e ° ° ° ° O 

0 O ee ° ° e ao) z. 20) O ° » e - ° O 
ry ~ O 


tee ) = 
k? BK | 
eee &, 6 Sy O_ - es > ee” tC) 


02 
+o 
03 
NO 
(A 
tt 
_ 
O 
O 
S) 


O Py e ° ° ° ) s O SB e . ° Bo 6) ©0| 


In the case of Equations (5-6) and (5-7), the constant zeta loci can 


be obtained from Equation (5-10) by making the following substitutions: 


a (-1) a, =) (5-11) 


tho 
2 cee, aU, CF ) (5-12) 


ja a 





The following example will demonstrate the technique: 


Example 2 


Consider again the symmetrical parallel-tee network snown in Figure 
2, Section 3. For the zeros of the transfer of this network, Equations 
(5-6) and (5-7) become: \ 
3 | 2 , Ee 
-XB WU, F) + 2AWUj (CF) - 2BUK FT) + U0 F) = 0 (5-13) 


-(f w) *u,( % y+ 26 wu 1% ) - BRE )+U,0¥)=0 G-14) 


Multiplying these equations ee uj); and 1 and forming the determinant 


204, 3 ), we “obtain: 


|xeU5 | +29U,  - 280, Uy 0 
0 a6, 429U,  -2 FU, ae 
O O eels veau ,, Teen 
2(X,B5¢) = iL v20U_, 280, - _, | 0 
0 -aBU,  +28U_, -26U_, ws 
1 ee 0 -4QUy  ¢2QU_,  -20U_, 





Substituting the appropriate functions of zeta for the UL. in Equation (5-15) 


and expanding, we obtain the equation for the constant zeta loci: 


~ 46 “(49 7-1)°[(-32 4 416) > + (Gh Fs16xF* - 3297 - 8K) B+ 
(5-16) 
(-64¢° 5° 4 ug? - 1267 97 + x*)] = 0 
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When the two equations (5-8) and (5-9) are of the same order so that 
m= on, Sylvester's (2n x 2n) determinant can be reduced to an (n x n) deter- 
minant by using Bezout's Method of Elimination (7). To demonstrate Bezout's 
method, consider Equations (5-8) and (5-9) with m = n = 3. This pro- 


duces the equations: 
3 


7. 
B(EODP2 fy 4 fi Met tm + £, =O (5-17) 
g(W)=2, W +g, WH +2, wi +e =0 (5-18) 
3 2 1 0 
Then: 
i5e@ - gf = 0 (5-19) 
(f£,0 + f,)8 - (8,4 +8,)f = 0 (5-20) 
2 2 : 
(f£,W° +f,W + f,)g - (g,W° +8, W +28,)f = 0 (5-21) 
can be written in matrix form: 
(£485) (£,8.) Ga oo 
382 38) 350 
(f£,8,) ((£485) a (f£,8,)) (£58) W) = (5-22) 
| (£38) (Eg) (£185) | oh 
Where (£585) denote (2 x 2) determinants: 
f. ay 
f.g.) = 5-23 
( 18) z, g, ( ) 


In order to eliminate ( from Equations (5-17) and (5-18) by Bezout's 
method it is necessary to force the determinant of the coefficient matrix 
in Equation (5-23) to zero. Thus the constant zeta loci for the case 


m =n = 3 are given in non-parametric form by: 


24 


—— 


**= el = 


a Gai 





(£385) 


Z,(f, 58) = 1 (£,38,) 


(£38) 


For cases where m 
developed by forming n 
such a procedure gives 


Zeca loci: 


(2.8) (f 
(2) S5) (2,6. ) 
ep et) ~ |. r 
(2.8, ) (2,85) 
\= £85) Ga 


(£,8,) ; (£485) 
((£,8,) + (£,8,)) (£85) | =O (5-24) 
(£,85) (£, 85) 


Hl 


nm pm 3, a determinant similar to (5-24) can be 


ee 
ees 


— 
4 


equations as in (5-19) - (5-21). Whenm=n=4 


the following determinantal equation for the constant 


ae, (2,8, ) (2.35) 
# (2365)) (2y89) + (258,)) (4589) 
+ (f56,)) ((2385) + (2567) (cecil esas 
5p! 235) (2465) 





The following example will demonstrate how Bezout's method is applied 


fo obtain the Gonstance 


Example Lil 


Let us again consider the symmetrical parallel-tee network. 


zeta loci: 


Equa- 


tions (5-13) and (5-14) yield: 


» 


“ABU, 
“AB Uv, 


59 


+-2 BU, 


+26U_, 


= -2BU, 


= -28U_, 


25 


aT 80 











Substituting the appropriate functions of zeta into these expressions and 


forming the (2 x 2) determinants for the (£585) we obtain: 


(£485) = +2% B47 -1) 

(£485) = -4% 6° F (4% *-1) 
(£48) = (4X A77?- 4 B (4S 7-1) 
(£,8,) = 48°F *-1) 

(£58) = -4 OF GE 2-15 

(£,8)) = +28 (43 “-1) 


When these are substituted into Equation (5-24) we obtain: 


+28 -485 Geet 
20%, 8,5) = KB 747 7-1)" -4% 8% (4% 77 +4B-%) -43 |= 0 
(4k FZ *-k) 45 +2 


Which when expanded gives the constant zeta loci in non-parametric form: 


Z5(4 Bf) = xP R71)? [(-32% 5 7416) 8? + (ue F416 % 72-3227 - 


Bx)B + (-6407 57 4 uBar et -120%7 3574 *)] = 0 
Observe that the determinant obtained by Bezout's method has the op- 
posite sign from that obtained by Sylvester's method. Sylvester's method 
produces the resultant of the functions f( @ ) and g( W& ), while Bezout's 
method produces the negative of the resultant (8). Since multiplying both 
sides of the determinantal equations by -1 does not introduce extraneous 
common factors of f( ™ ) and g( WJ), this sign difference can be ignored. 


It would appear at first that Bezout's method would be easier to apply 
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ou 


than Sylvester's method because of the 


mMinancs thet it affords. 


Peewection in the order of the deter= 


However, Sylvester's determinant has many null 


elements as well as functionally simpler elements than Bezout's determinant. 


It is the author's opinion that Sylvester's method is somewhat easier to 


apply. 


Equations of the Constant Omega Loci 


Using Theorem (II for the cases j = 0, j = 1, -, J = (n-l) we 
eptain: 
. } Ic 
c 
"=O 
k k 
ore ee = 
K=O 3 
(5-26) 
“") 
k k: 
» (-1)"a, © U(F)1) = 0 
kK=0 
These equations may be written: 
0 0 1 iL n n - 
(-1) a, w& Uy + (-1) ay aU, - + (-1) a De = 0 
O 0 1 oh n n 
= iw - 7 - =i) — - 
(-1) ay Uy + (-1) a, W Uo + + (-1) aw U1 0 (5-27) 
Q° 0 1 n ig 5 
(-1) ay w Yf-i" (-1) ay +(-1) aw U, = 0 


wus. iy ae 


Making use of the relation UC? ) = - U_¢ i) and Uy 


the set of equations (5-27) we obtain: 


ee 


i (-2)"a,w" ee 
ee 5) (a) eos 
3 Wa w'a(-1)" "a 4 vo 

pay ia _# en 


= Q and rearranging 





4 fod 
Cee ee ee 
n 1 
O U, 
=O (5-28) 
fe U 
S 
O wo Thee 
(-1) anid. | | = 
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ne? eae aaa 


The set of equations (5-28) are linear and homogeneous in the quantities 


U Thus, for a non-trivial solution to exist, the determinant of the 


I” 
coefficient matrix in (5-28) must vanish. This condition specifies the 


relation nacessary among the a, anda) which specifies the constant omega 


kk 


Leroi 
- i 
(-1) a wt ee (-1)"a wo” 
ol n 
2 
W(a,,@) = ie 
on on n-2 n-2 (5-29) 
((-2) aw ==) RD eae O 
n=1 = 36 
(-1) Gp eee (-1) Bobo 


Example IV 


Let us apply this procedure to obtain the constant omega locus for 
the zeros of the symmetrical parallel-tee network. From Equation (3-2) 


ERe Zeros Ofstnewtrenster function for this network are Ene roots or: 


wBp + 2@p + 2Ap+1=0 (5-30) 


Applying theorem II, we may write as in Equation (5-26): 


eat -'% Z 2 é) 3.) ae 
G1) (26) du, | (-1)*(28) WU, + (1) (a8 ) OU, 0 
0 G Z 
(-1)°2) w UL, + (-1*(28) wt, + (-1)°( xB) wu, = 0 (5-31) 
OL 0 ae 1 1 a 3 5 2 
(1) wu + C1 (2p) wu, + CDMA) au, = 0 . 
Using the relation U ae = - Uy. these equations can be written: 


28 





(-1)'(28 Ciece 0 (-1)°(@@ Ww? 


(-1)7(28) 07-(-1)"°C) w} (-1)2( 8 yw? 0 


[(-1)° (xP )w?-(-1)'(28)o']} —-(-1)"(1) 0 0 





Setting the determinant of the coefficient matrix to zero we obtain: 


- 2A +2 Bu” - ot ay 
WA, 8,0) = (2 w*-1) -& Bw? 0 =0 (5-32) 
(28W - x BW) -1 0 

Which, when expanded becomes: 
WOK ,8,W) = KA(a7 Aw - 2K —e7wr + 2BW*%-1)=0 (5-33) 
Thus the equation 

(or 74° - 2 XW") : a (+2 0) - 10 (5-34) 
gives the constant omega loci except where & = 0, & = 0 or W= 0. 
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6. Equations of the Constant Bandwidth Curves 

In this section, equations will be developed for displaying curves of 
constant bandwidth upon the parameter plane. 

Let us define what is meant by a constant bandwidth curve: 

"A constant bandwidth curve for G(j Ww) = M is a curve drawn upon the para- 
meter plane which specifies the relation between the parameters necessary 

if the transfer function G(s), which is a function of the parameters, is 

to have magnitude M at the real frequency 4) b' 

Once such curves are obtained for selected values of M, the frequency 
response G(j o) ) may be sketched. Alternately the constant bandwidth curve 
corresponding to some M and W specification can be drawn, and from this 
curve values of the parameters can be found which will guarantee that the 


specification is met. 


Consider the following rational transfer function G(s): 


s" “ip 3” +t ese teas ad 
OS ol Se a ee (6-1) 


rs + fr S + wee + 
a r,5 + To 


Where the coefficients of Q(s) and R(s) are of the form: 


q 


ex +f + +h. 3; p.™m0, 1, ..% 6-2 
A 4 ug 8,48 aes? me (6-2) 


ul 
© 
_ 
hd 
e 
= 


ry) ay + b 6 + c. XB ~ qd. ; k (6-3) 


Thus, G(s) may be written: 


m 
Va eA+ f A+ X +h sP 
‘ P “ph ep (s p? 


G(s) = (6-4) 


J 


7 k 
(ag+b A + cM +d) s 


ll 
© 
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Consider now only real frequency, s = j 4 , this can be written: 


G(jW) = 

= kp op — k(p-1)  p 

Lu; (-1)" wi (ed +£ p+8,x@ th.) + 3 dy (-1) uw“ (e,4 +£,6 +8,«6 th,) 
even odd 


=. n 
my (-1)* w* (a,x +b, 8 te,x td.) + 5 Bs (18D ya, +b, B+c,~6 +d.) 


even odd 
(6-5) 
Now let W= W,, some fixed value of &, and also define: 
n 
A = > 7 we peetc.afom Be. C.D (6-6) 
r b k Yr r re 
k=0 
even 
n 
A = >: (1) oe ; etce forss,, C.D. (6-7) 
i bk i il i 
k=1 
odd 
m 
E = a (ye? ae e ; weve... for F ,G “a (6-8) 
r = b p Gor er 
p=0 
even 
m 
ES = ey Gye re s/s etc. for F,, G., H. (6-9) 
ol baa iL yee | 
odd 


Using definitions (6-6) thru (6-9), Equation (6-5) can be written: 


(ME + BF So (3 eS te H.) + jCOE, + OF, +A G; ate H,) 
G( 3 wi) = (of A+ QB. + #3 C. s. D.) + jCXA, a2 (6 B. eee C. =e D,) (6-10) 


@ 


Using (6-10) and defining the magnitude of G(j wW) as M we obtain: 


oul 





2 2 


2 2 
M= |G(jM,)| ~ = 2 2 
b (x A+R B+ of 8 C+D.) +COCA QB + x /3 C+D.) 


Which when expanded becomes: 
M = Q(x 46) (6-12) 
Where; 
Me .2).).2 2 2 2 
aa, pmo gp (G+ Gi) + 2k BEG +£G,) + 2%8 (FG +F,G,) + 
2, m2 2 2p 2 2 
XX (E “+ E;) + Bp am + F.) + 2x6 (EF + EF.) + (6-13) 


' ; /) eee?) 
2X (EH + EH) + 26 (FH + FW.) + (+ Hy) 


2105 > ane | 2 
O(X, 8 ) Sica (C 4. C.) + 2X8 (Ac. . A.C.) TP ates fe (BC. + B,C.) + 


ly ae ae ee Z 
of (A. A.) - @ (BL + Bs) + 2X8 (A BL + A.B.) + 


) 2 2 
2x (AD. + A.D) + 26 (BD. + B.D) + (D. + D;) (6-14) 


Equation (6-12) can be written: 


Dog) - (MO (X,@)) = 0 (6-15) 


Substituting (6-13) and (6-14) into (6-15) and expressing the resulting 


equation in the form of a quadratic in B we obtain: 
(MR,-R,) @ 7+ QW, - WB + OPV,-V,) = 0 (6-16) 
- 2 1 a m2 


Where: 


a ee 2 22 2 
R, = (B + B.) +2 (BC. “ B,C,) + o (Cc. a C.) 
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-*) 
i! 


2 A Loe: 2 
cr + FD + 2“ (FG. + FG.) + OX (Gc. + G.) 


= 
i 


2 
2(B.D_ + B.D.) + 2X (AB +A,B,+C D4C.D,) + 20° (AC +4A,C,) 


2(F H + FH.) + 2X(E F +E.F.4G H+G.H.) + 2%(E G +E.G.) 
rr Lc 2 Mesias Ul ) gee oy iL r rx Yr YL Ler pe li 


=x 
ul 


2 2 2,.2 2 
(D- ~ D;) + 2X (A D_ + A.D) +X (A, -- A.) 


< 
| 


< 
" 


2 2 Lee 2 
of 
(HT + H.) + 2 (EH + EH) + x (E + E;) 


We may now solve the quadratic equation (6-16) for p as a function of &, 


Wy» and M. This solution is: 


2 2 
-(W, -W, z (MW, -W,) - 4(M°R, -R,) (MV, -V,) Gis 


Pigs * 2 
: 2(M R,-R,) 
The constant bandwidth curve corresponding to some desired M = My and 
W, = WwW can be computed from (6-17) by a digital computer. The 


b bO 


procedure is to fix M and &, in (6-17) at the desired values, increment 


b 
over the range of interest, and at each value of & to compute the corres- 
ponding values of  . The two roots of (6-17) give the two branches of 

the desired constant bandwidth curve. The reader is referred to section 10 


where constant bandwidth curves are computed and discussed for a loaded and 


null adjusted parallel-tee network. 


33 








‘as Equations of the Constant Q Loci for Second Order Systems 
Recent interest in the Q of second order RC networks, accompanied by 
difficulties in applying the conventional definitions of Q, led Morris (4) 
to define Q mathematically as a parameter of the characteristic equation. 
He rearranged the second order characteristic equation into the form (7-1) 
which defines Q: 
Ze 


Ts +—yst1=0 (7-1) 


Dutta-Roy (9) has pointed out that Q so defined is equivalent to an 
earlier definition proposed by Bolle (10). Bolle's definition for the 


second order characteristic equation: 


2 
as + a,s + ay = 0 Tae 
si: 
(aay)? 
: ea (7-3) 
1 


In this section equations will be obtained which will permit loci of 
constant Q to be displayed upon the parameter plane for variable parameter 
second order characteristic equations of the form (7-2). It is assumed 
that the variable parameters appear in the coefficients of (7-2) in the 


following manner: 


b,x + cf ci de a eo (7-4) 


“k 
k = 0,1, and 2 
Squaring both sides of Equation (7-3) and rearranging we obtain: 
Q7a - a.a, = 0 (7-5) 
1 Oe 
Substituting (7-4) into (7-5) and writing the resulting equation in the 
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form of a quadratic in , we obtain: 
2 2 2 Zz 22 2 
[(Q d,-dod.) + (2Q cd, -cod,-cod aX + (Q C1-fpc,) ] (3 + 
[(297b d.-b.d,-b,d 9 : + (297b Cc +2074 f. -bee “bec” -d fd T-) o + 
tr GZ 2°70 11 1livwz0 C2 ru 2 


(29°c, £,-e4£.-epf,)! 3 “ 
22 2, ,.0°b f tbaf bE + (Q° FFE} =O (eb) 
[(Q=b]-byb,) & + (2) =r Oro 2 1°02 
The desired constant Q loci can be obtained from (7-5) by use of a 
digital computer. The procedure is as follows: 
ee Fix Q in (7-5) at a value for which the locus is desired. 
2. Increment &% over the range of interest. 
oS: For each value of & , solve (7-5) by the quadratic formula. 
4, The two solutions of (7-5), B, and rm) »? plotted as ordinates 
on the parameter plane with &% as abscissa define the two 
branches of the desired constant Q locus. 
De Repeat steps 1-4 for each value of Q for which a locus of 
constant Q is desired. 
For examples of this procedure the reader is referred to Section 10 


where several loci of constant Q are obtained for a variable parameter, 


loaded, and null-adjusted parallel-tee network. 
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obs Sensitivity Analysis by Parameter Plane Techniques. 
a. Introduction 


The polynomial equation 


>on a, = 0 (8-1) 


where : a = £(q,»d5>-+->4,) 

has n roots. Each of these roots is a function of the parameters Gy? Io» 
oe: eee ae It is the objective of this section to show how each of 
these roots vary as the parameters qd change. 


Let us denote the ae complex root pair of (8-1) by: 


r 2-5, 8, ti Wr F’ (8-2) 


and the ;= real root of (8-1) by: 
re = OO .. (8-3) 


Let us now define two types of sensitivity; macroscopic sensitivity 
and microscopic sensitivity. 

Macroscopic sensitivity is the sensitivity of the roots (8-2) and (8-3) 
to large changes in the parameters q-° These sensitivities may be defined 


by: 
ec IAN ? C4 
if mK > ( 9 =) 4; =0, ifr 


S 


4), 
el SIR i)! S / > 429; =0, i#r (8-4) 
sie (AT DM AG +) 
j,° j 7 Lg aH =0, ifr 
sf denotes the change in the damping ratio 3 i of the complex root pair 


a) . 
r. due to a finite change in the parameter q-° S 7 denotes the change in 
9 
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the undamped natural frequency of the complex root pair r. due to a finite 


G? 

change in the parameter qd Ss. denotes the change in the real root 
9 

‘4 due to a finite change in the parameter qd 


Microscopic sensitivity is the sensitivity of the roots (8-2) and (8-3) 


to infinitesmal changes in the parameters qd These may be defined as: 


Meee Py! OF 

Ld,” _= 3 

. dw ;/ dg r (8-5) 
e 

Sir * ao~ 5199 13 


The interpretation of these microscopic see CANES except that the para- 
meter changes are infinitesmal rather than finite. 
love Obtaining Macroscopic. Sensitivity from the Parameter Plane 
The procedure for obtaining macroscopic sensitivity directly from 


the parameter plane is best demonstrated by an example. 


Example V 


Let us consider some third order polynomial equation: 


3 2 = - 
248 + a,s tas ta, = 0; a, = £,.00%5°B) 


Let us also imagine that we have obtained the parameter plane for this 
polynomial and that it appears as shown in Figure 3: Let us presume that 
we have selected “= & _ = 3.8 and = 8 9 = 2-4 in order to place a pair 


0 
of complex roots at Zz = 0.2 and W = 0.7 and that this choice results in a 


real root at 06 = - 3.0. We wish to determine the macroscopic sensitivities 
° z WJ w) o o , —_ 
Sx A Se ; :” ; Se : Sy » and Se » for an increasé of one unit in 
both & and f@ . From Figure 3 we determine: 
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Fig. 3 Macroscopic Sensitivity from the Parameter Plane. 
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More important perhaps, than these numerical values of the macro- 
scopic sensitivities is the qualitative, but enlightening information 
about the sensitivity of the roots which can be obtained by careful in- 
spection of the parameter plane. Bethe constant zeta, omega, and sigma 
loci have been drawn for constant increments in each of these quantities 
respectively, then one can immediately see that the most sensitive areas 
are those in which the loci are tightly spaced. This is like looking at 
a contour map of terrain- steep areas are highly sensitive areas. 

come Obtaining Microscopic Sensitivity from the Parameter Plane 

Microscopic sensitivity can be obtained either by the method of 
Kokotovic and Siljak (5) or by an alternate method to be developed here. 
Kokotovic and Siljak derive relations for sensitivity from the basic para- 
meter plane equations. The new method, to be described here, obtains 
sensitivity relations directly from the non-parametric equations of the 
constant zeta and omega loci which were developed in section 5. In those 
cases where the equations of these loci have already been obtained by pre- 
vious calculation, sensitivity relations can be more rapidly developed from 
these loci equations than by the method of Kokotovic and Siljak. As back- 
ground, the method of Kokotovic and Siljak will be discussed first, and 
then the method of obtaining sensitivity relations from the non-parametric 


equations of the constant zeta and omega loci will be introduced. 


(1) Method of Kokotovic and Sil jak 


The basic parameter plane equations are: 


f= > (-1)ta, od “U,_1( ¥) = 0 (8-6) 
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n 


k k 
g = 2, (-1)"a, wU( J) = 0 (8-7) 


Functionally these equations can be written: 
f£(F ,W, a) = 0 (8-8) 
g( 3 ,W ) a, ) = 0 (8-9) 
Consider that a and @ in (8-8) and (8-9) are implicit functions of the 
a,. The a, are in turn explicit functions of the qd 


k 
Thus: 


9£/9q, = (9 £/0W,)( IW ,/9 4,)4+(9 £/ OF) (OF /2 a,)+(9 £/ a) (0 a,/2q,) 


(8-10) 
dg/dq, = (Og/dw,)(dwW /O a,)+( 0 8/ O% )(IF19q,)+(A8/0 a,)(da,/dq,) 
(8-11) 
But from (8-6) and (8-7): 
n 
df/0W = (1/ 4.) Ginna k ae 
rt i 2 ay a)! U,_76 Z,) = AL (8-12) 
n 
» = Fe rae 
D1 be, = » (-1) "ay eu mie, A, (813) 
n 
Oe/ dm, = (/ ,) a (-1)"ka, wu, ¥.) 2 Bt (8-14) 
n 
dg/O¥, = a (-1)"a, WU) 2 By (8-15) 
k=0 ‘ 


Substituting (8-12) thru (8-15) into (8-10) and (8-11) we obtain: 


40 





n 


> £/q, = (/a .)¢ dw jd q_) 2, (-1)* tka, w/ Uy + 
n 
(OF 104.) = (1 ey ex + 
k=0 . 
n 
ba (Da /Oq)(-1)* we, 4 = 0 <eare? 
n 
dg/dq, = (1s .)¢ dw Ue q.) 2 (-1) ha, w)-U,, + 
n 


(OF 04.) 2 (-1)a, wus + 


n 
2 (Ja /dq_)(-1)* w.“U,, = 0 (8-17) 
Let: 
n 
&. = e (Ja/dq_) CDT yD (8-19) 
n 
Ce = pa (Jaq) Gaye Ww, (8-20) 


We may now rewrite (8-16) and (8-17) as: 


(Ai) @ ,) JW./dq, + (Bi) OF /0q, = -c} | (8-21) 
i i i 
G5!) day Jod. + G5) Oe vem a -coe (8-22) 


Equations (8-21) and (8-22) are simultaneous equations, from which the 


desired sensitivities ( OF Je qd.) and (dW JO qd can be obtained. 
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d. Sensitivity Relations from the Constant Zeta and Omega Loci 

In Section 5 it has been shown how the equations of the constant zeta 
loci 2a,» ys ) = O and the constant omega loci W(a,» 4a) ;) = 0 can be 
obtained in non-parametric form. From these equations we can obtain the 


relations for the microscopic sensitivity directly: 
Z(a,; x: ) = 0 (8-23) 


W(a,, w ) = 0 (8-24) 
from which: 


OF Jda, = -(Dz/da)/(0 2/ 98 ,) (8-25) 


providing dz/ 2) x i # 0 
dw 7d a = -(9 w/a a)/ (2 W/ Ju.) (8-26) 
providing dW/d w i # 0 
For an example of sensitivity computed from the non-parametric equa- 
tions of the constant zeta and omega loci the reader is referred to 


Section 11 where root sensitivity to parameter variations in a lattice net- 


work is discussed in detail. 


42 








2. Parameter Plane Analysis and Design of the Symmetrical Parallel- 
Tee Network 
a, Introduction 

The symmetrical parallel-tee network shown in Figure 4 has 
been the subject of several papers (9,10,11). The reason for this interest 
stems primarily from the fact that it is possible to obtain complex zeros 
of transmission using the network. These may be placed anywhere in the 
left half s-plane and as far as 30 degrees from the imaginary axis into 
the right half s-plane by suitably choosing the network elements. From 
the transfer function of the network, Equation (9-1), it can be seen that 


the poles and zeros of transmission are functions of the parameters of and 


B 


Eo(P) = Bp +28 p+ 2B p+] (9-1) 


Pi) Bp + (K42B + 2B )p (428 +2)P+ 1 





where: p = RCs 

Barker and Rosenstein (12), using root locus techniques, investigated 
this network. Their analysis, while complete and well done, was compli- 
cated considerably by the inherent single parameter capability of the 
root locus technique. This section approaches the same problem with a 
two parameter method of analysis and design - the parameter plane. 

By employing parameter plane techniques, graphs are prepared for 
the polynomials which form the numerator and denominator of the transfer 
function. These graphs, with @& as abscissa and B as ordinate, display 
directly the roots of the polynomial which result for any choice of the 
parameters & and iE . From these parameter planes for the network, the 


regions of the s-plane in which poles and zeros of transmission are possible 
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are immediately evident. Likewise values of the parameters & and g8 
necessary to place either poles or zeros at some desired point within these 
regions can be directly read from the curves. The parameter plane is there- 
fore a tool of both analysis and design. 

Loci of constant bandwidth are drawn upon a parameter plane for this 
network. These consist of a family of curves along any one of which the 
frequency at which the response of the network is down N db is constant. 

For any choice of & and a the frequencies at which the response is down 
N db can be read directly from these loci. Loci for the conventional -3db 
bandwidth are presented in this section. 

Thus it is possible to select values of the parameters o& and I 
which guarantee suitable s-plane locations for certain poles and Zeros 
while indicating the bandwidth determined by this choice. Conversely & 
and a may be chosen to guarantee upper and lower -3db frequencies while 
offering a selection of poles and zeros which may be read from the curves. 

The curves which are presented in this section are frequency normalized 
and hence universal. They may be used to investigate any symmetrical paral- 
lel-tee network of the form shown in Figure 4. 

The advantage of the parameter plane over other methods which have 
been used to analyze and design this network are the rapidity and ease with 
which the capabilities and characteristics can be determined from the curves 
and with which the parameters may be selected in order to produce desired 
network behavior. | 

b. Approach to the Problem 
The objective of the analysis to follow is to determine 
in what manner the poles, zeros, and -3db frequencies of (9-1) change as the 


parameters o& and are varied. This can be accomplished by preparing three 


A4 





parameter planes; one each for the poles, zeros and -3db bandwidth of (9-1). 
Since some of the equations of these parameter planes have already been 
developed in earlier sections, they will be recalled as appropriate. 
S Poles of the Transfer Function 
Since the network shown in Figure 4 is passive RC, all of 
the poles will lie along the negative real axis in the p-plane. Let & 
denote the value of p at which such a pole of (9-1) occurs. We may thus 


assert that poles occur at those values of @ which satisfy (9-2): 


MOT + (2K 4+ 2AB th)E* + (ZB+ 2K 4+ 2TH+1=0 (9-2) 
Rewriting (9-2) in the form of a plane algebraic curve in the &%- 
plane we obtain: 

(e242 vr“) AG + (26° +27) t$ (e740 )e + (20+ 1) =0 

(9-3) 
Equation (9-3) represents a family of hyperbolas in the&X- 8 Sian One 
hyperbola occurs for each selected value of GW. Since (9-2) is satis- 
fied by three values of @&@ for each choice of A&A and 8 , it follows that 
each point in the parameter plane will have three hyperbolic branches pass- 
ing through it. 

The curves represented by (9-3) can be hand calculated for selected 
values of @ , but since this is a tedious process, a digital — was 
used to obtain them. The resulting curves which specify the pole values 
for any choice of the parameters are shown in Figure 5. One enters this 
figure with values of & and é and then reads the values of the three 
poles directly. Alternately one may select values of & and 6 to produce 


some selected pole values. 
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di. Zeros of the Transfer Function 
In this instance the problem is much more interesting since 


complex zeros are possible. The zeros are the roots of: 


3 
of (3 P +2 Bp +26 p+1=0 (9-4) 
Since and /3 are presumed real it can be asserted that (9-4) has at 
least one real root. As before let p =O in (9-4) and rewrite the re- 


sulting equation as a plane algebraic curve in the & - S plane; 


o3ag t(2o° +20 dB +1=0 (9-5) 


Equation (9-5) also represents a family of hyperbolas in the & -f plane. 
These are shown in Figure 6. Note that in some regions it is possible to 
have three real zeros. 

In describing the location of complex zeros of (9-4) the zeta and omega 


notation familiar to control systems theory will be used: 


s = - FW + jw yl- 3? (9-6) 


The non-parametric equations for the constant zeta and omega loci which were 
developed for this network in Section 5 will be employed. The equations of 
these loci are repeated below. For their derivation the reader is referred 


to Section 5. 
UX, BF) =- KB F* -1)? ((-324K 77 + 16) B+ (OAK + 


16 & 7 -32 ie 84 )6 + (-640 77 © 4 “gay * - 12K*x*? +a 4) = 0 


(9-7) 


18 ay? («Qu - 2 Ba)" + 2 Ba? - 1) =0 
(9-8) 


WOH , Bw) 
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A digital computer, programmed as described in Section 5, was used to 


obtain the constant zeta and omega loci which are shown in Figure 6. 


e. -3db Constant Bandwidth Curves 
One may obtain the desired -3db constant bandwidth curve by 
employing the procedures which were developed in Section 6. Referring to 
that section we note that several quantities must be substituted into 
Equation (6-17) in order to obtain the curve which we desire. In order to 
facilitate substitution into (6-17) the following quantities are tabulated: 


From the transfer function (9-1) using the notation of Equations (6-2) 


and (6-3): 
7; 5° ae Se =e Ss SS SE 
0 0 0 0 IL 0 0 0 1 
1 jO, 2 0 j20, 0 j20, 0 0 
2 2 2 2 
2 ed - 2) - 2H). 0 0 - 2W) 0 0 (9-9) 
3 0 0 - ja), 0 0 0 - jw 0 


From Equations (6-6) thru (6-9) using values from (9-9) one obtains: 


A B c D E F G H 
2 
Real -Wi, ~20, -20* 1 0 - 2u 0 1 
(9-10) 
3 3 
a) = 
Imag @)y, 2u) 0), 20) 0 2u) wr 0 
From Equation (6-16) using the values from (9-10) one obtains: 
4 Pe 4 Zane 2 6 
= 40) W HK w) 
a vee 5 + 4 , rae Ota w, 
i 2 2en0 4 
= 2) 
R, aa, + _ on D4 Wy 
2,4 2 
Ws = 2 x Ww ~ 4 Ww) (9-11) 
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2 


W, - “40, 

2.4 Ze 2 2 
sabes Oi t+ x MW, + 2XW) + WwW 
Vo = 1 


The quantities (9-11) are used to digitally compute values of Equation (6-17) 


corresponding to selected values of &, with M = 1/ J 2 as required for the 


b 


-3db curves. The resulting curves are displayed in Figure 7. 


f. Interpretation of the Parameter Planes 


The parameter plane for the poles of the transfer function 
which is shown in Figure 5 reveals that for 0 € &, 4 < 2.5 the three 


poles are confined to regions along the negative real axis as shown below: 





The parameter plane shown in Figure 6 for the zeros of the transfer 


function reveals the following information regarding the zeros: 

(1) The positive * axis coincides with the 5 = -% locus 
indicating that zeros may be placed as far as 30 degrees from the imaginary 
axis into the right half p-plane. 

(2) The zeta = 1 locus separates the region of 3 real zeros 
from the region of 2 complex conjugate zeros and 1 real zero. 

(3) The zeta = 0 locus is a straight line with equation 
A =4( . This indicates that if one is interested only in those cases 
where a pair of conjugate zeros are placed on the j@ axis, then he may 


substitute this linear relation between the parameters into the network 
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eguations thus reducing the problem of analysis and synthesis to one of a 
single parameter. 

(4) The zeta = % locus coincides with the omega = 1 locus 
indicating that an additional degree of freedom exists which may be used 
to specify other poles or zeros in addition to the specification that a 


pair of complex conjugate zeros occur at y= 5, G= 1. 
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Fig. 4 Symmetrical Parallel-Tee Network 
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Fig. 5 Poles of the Symmetrical Parallel-Tee Network 
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Fig. 6 Zeros of the Symmetrical Parallel-Tee Network 
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Fig. 7 -3db Constant Bandwidth Curves for the Symmetrical Parallel- 
Tee Network 
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10. Loaded and Null Adjusted Symmetrical Parallel-Tee Network - 
Parameter Plane Analysis 
a. Introduction 

The null - adjusted symmetrical parallel-tee network provides 
a means of obtaining frequency selectivity without the use of inductors. 
Hence it is a very useful network at low frequency where inductor size be- 
comes unacceptably large as well as in integrated circuits where the current 
technology does not permit inductors of reasonable size to be integrated. 

A number of articles have considered the zero source - infinite load 
resistance case (9,13,14), and a few have considered the effects of other 
finite combinations (10,11,12). These latter articles generally consider 
either the source or resistance to be fixed and the load resistance to be 
variable or vice-versa: 

This section considers the problem from a different point of view and 
with a two parameter tool of analysis and design - the parameter plane. 

The two parameter nature of the problem is evident by inspection of Figure 


8 and Equation (10-1): 


2 


ace) _ Eee (10-1) 
Ei(P) (4B +l +B) p24(4 HB 44% 44642)p + (K +h+ 2) 





where: p = RCs, A= R_/R , f= R,/R 

Parameter planes will be obtained for this transfer function. They 
will show the poles, Q, and frequency response characteristics associated 
with this network which result for any choice of the parameters & andf@ with- 
in the range 0¢ ~,A@K< 3. Alternately 4 andf& may be chosen from the 


curves to obtain desired characteristics within the capabilities of the 
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network. 

Loci of constant Q discussed in Section 7 are also obtained. These 
loci show directly what value of Q results for any selection of parameters 
o and B® . Conditions are specified which permit Q to be related to the 
bandwidth of the network. 

Constant bandwidth curves are obtained as described in Section 6. 

For any choice of parameters these curves immediately give the upper and 
lower frequencies where the response of the network is down N db from the 
response at infinite frequency. Curves are presented for several values 
of N, thus permitting the frequency response to be sketched for any chosen 
values of source and load resistance. 
b. Poles and Zeros of the Transfer Function 

The parameter plane graphically displays the roots of a 
polynomial whose coefficients are functions of two variable parameters. 
In the general case these polynomials (which may be the numerator or denom- 
inator of some network function) can have complex as well as real roots. 
In the case of Equation (10-1), however, the following observations are 
made: 

The zeros are fixed at me =- 1] 

The poles are all real and lie on the negative real axis 

in the p-plane since the network is passive RC. 
Thus it is necessary to investigate only the poles of (10-1) and this is 
considerably simplified by the fact that the poles are all real. To obtain 
the parameter plane for these poles let p =Q bea pole of (10-1). Making 
this substitution into the denominator of (10-1) and rearranging the result- 


ing equation into the form of a conic section in the &-f plane one obtains: 
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467440) + ema 4: 1)A + (¢* + 474198 + (20+ 2) = 0 
(10-2) 
This equation represents a family of hyperbolas; one hyperbola for each 
value of © . These are plotted for several values of O and are shown in 
Figure 9. One may read the poles, which result for any choice of & and, 


6 , directly from the curves. 


Example VI 


Assume that a parallel-tee notch network has been designed to reject 
60 cps. The notch filter has been designed so that R = 2660 ohms and C = 1 
microfarad. The network is driven from a source resistance of 5000 ohms 
and drives a load resistance of 2500 ohms. It is required that the poles 
of the transfer function be obtained. From Figure 8 one obtains: 
oe = 5000 = 1.88 
2660 


B = 2500 = 0.94 
2660 


From Figure 9 the normalized poles are read: 
gq, = -1.78 
Or = -0.27 


Which, when un-normalized become: 


Ts ee = Be = -670 


2.66 x 107° 


-0.27 


on 3 


———=— = -102 
2 2.66 x 10 


The network transfer function is thus: 
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2 
E(s) = (3 [ - + wee ] 


E,(s) 4adf +X+ 8 (s +o) Ge y) 
= 0.95 [ a ees ; 


(s+670) (s+102) 


The curves of Figure 9 can be used to place one of the poles at a de- 
sired location while maintaining some control over the location of the other 
pole. This is accomplished by changing %, #, or R. Changing R and 
assuming that C is also changed appropriately to keep the zeros fixed changes 


both & and f. 


Example VII 


Suppose that in the previous example one is able to change the source 


resistance and desires to do so in such a manner as to cause the pole at 


s = -670 to move to s = -750. Thus the new normalized pole location be- 
comes: p, = -2.00. Since the load resistance is fixed at 2500 ohms, 
remains 0.94. From Figure 9 note that the G = -2.00 locus intersects the 


6 = 0.94 coordinate at & = 1.07, and that the other pole is thus located 
at py = -0.335. When un-normalized these pole locations become; 


sy = -750 


Ss -126 


Z 


The new transfer function becomes: 


Es @r156 [ se + 3777] 
E.(s) (s + 750)(s + 126) 


The new source resistance required becomes: 


R, = & R = 1.07(2660) = 2840 ohms 


56 








Ce. Constant Bandwidth Curves 
By making the substitution p = j W as discussed in Section 
6, the real frequency transfer function is obtained for the loaded and null 


adjusted symmetrical parallel-tee network: 


E(w ) 


mor OKC AK GB) (CCX. ae) (10-3) 
E,C@ ) 


Where: 
3 


KC“ ,8) = 44 +x+ 7 


(10-4) 


2 
eee EE 


[(X +B +2)-(4K 8 +048) W°] til (408 Hoey 446 +2)] 


By inspection G(e& , f > 4) )—~1.0 as W- eam . At any selected frequency 


wa) b? and attenuation level M at that frequency, the relation 


G(x, 8 ,W,)=M (10-6) 


defines a constant bandwidth curve in the &-/@ plane. This curve specifies 
the source-load resistance relation requires for (10-6) to be satisfied. 
The attenuation level M can be thought of as the attenuation from the in- 
finite frequency value of E 6 2 )/E,(& ),; which is K(X, @ ): 
K(&%, ®& ) = Be ae 
E,(@ ) (10-7) 


In order to obtain the equation of the constant bandwidth curves Equation 
(10-5) is substituted into Equation (10-6) and both sides of the resulting 
equation is squared. This squared equation can be solved for 8 as a 


function of Q , Ww» and M: 
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-(MW, Wy) + [QMCW, -W)°-4 (MR, -R,) (HV, -V,) 2 ae 


2 
2(M R,-R,) 
Where in the notation of Section 6: 


Ridelwi 4+ 140 1 + GoiCetany 6005) 4 lo a-car4 wa 


ih 


we be + 1 Pde Gow ue) Slee Cw 2a 1) 


w 
it 


2 2 4 


WD Sie teceeemtdo” ail ameter eae colo) 


= 
il 


2X(4)t-2 W- +1)+ Beda or oe + 1) 


= 
il 


a oe + 4 (307 + 1) +x 2(w* + 14 W? + 1) 


< 
ul 


Vv, ™ x 7 w'-20" + 1) 


Assisted by a digital computer, Equation (10-8) was used to obtain the 
constant bandwidth curves for values of M corresponding to -3db, -6db, and 
-20db. These curves are shown in Figures 10, 11, and 12 respectively. 
These curves may be used to sketch the frequency response as indicated in 


the following example. 


Example VIII 


For the conditions specified in Example VI, namely & = 1.88 and 6 = 


0.94, the curves of Figures 10, 11, 12 are used to obtain: 


Db down from infin- Wr a) 

e . b 

ite frequency gain _l 2 
-3 0.50 2.65 
-6 0.61 1.84 
~20 0.901 gO bes) 
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The infinite frequency gain is: 
K(1.88, 0.94) = 0.095 
The zero frequency (DC) is: 
a = = 0.195 
E, (0) X+6+ 2 
The attenuation at W = 1.0 is infinite since the network of Example VI 
was designed for this condition. 
The above information may be used to sketch the frequency response 
of the network. 
ale Constant Loci 
As pointed out in Section 7, the Q of a second order network 
may be defined as: 
at geometric mean of oth and ae degree terms 
- degree term 
where these terms refer to coefficients in the denominator or the network 
transfer function. For the loaded and null-adjusted symmetrical parallel- 


tee network, Q becomes: 


one [(4¢8 4+6%+8)(¥4+6 + 2) 1? 


4 48 + 4% 44842 (10-9) 
Squaring both sides of Equation (10-9) and rearranging the resulting equa- 


tion into the form of a quadratic in B obtains: 


[16976 7 +(3207-4)e 4+(1697-1)) 8 7 + 
[(3297-4) %° si (4897-10) & a2 (1697-2) ] 8 7 


[(eg--1) > 4° 1eO--2)e Ph Go-) =O (10-10) 
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Equation (10-10) was used to obtain the constant Q loci displayed in Figure 
[Be 

Q defined in this manner is a mathematical parameter of the network as 
discussed in Section 7. It is not in general directly related to the 
frequency response, however, if the infinite frequency response is made 


equal to the DC response as explained later, then the following relationship 


holds: 
tO 
Q = i= oe (10-11) 
2 1 
Where: fo is the notch frequency 
fy is the lower -3db frequency 
f,, is the upper -3db frequency 


It has been shown (10) that the null adjusted symmetrical parallel- 
tee network under the conditions R. = Q, Ry = has a Q equal to 0.25. 
Figure 13 shows that Q may be made larger than 0.25 by adding source re- 
sistance. The maximum attainable Q is Ge and occurs at Re = Rp = 0./707R. 

e. Locus of Symmetrical Frequency Response 
Also displayed on Figure 13 is the locus of symmetrical 

frequency response. This locus defines the source - load condition re- 
quired for the DC response to equal the infinite frequency response. 
Cowles (10) has shown that this requires RoR, eS Ral ewee equivalently 
x 8 =i, Figure 13 shows that this locus also passes thru the point 
of maximum Q. Thus a loaded and null adjusted syatietrical parallel-tee 
network whose source and load resistances are chosen so as to maximize 
Q will also have a symmetrical frequency response. As pointed out above 
in sub-section d, the network under these conditions also is the most 
selective by virtue of its optimum Q. 
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Fig. 9 Poles of the Loaded and Null Adjusted Symmetrical Parallel- 
Tee Network . 
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-3db Constant Bandwidth Curves for the Loaded and Null 
Adjusted Symmetrical Parallel-Tee Network 
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Fig. 11 -6db Constant Bandwidth Curves for the Loaded and Null 
. Adjusted Symmetrical Parallel-Tee Network 


64 





{ 





2.8 


2.4" 


Fig. 12 -20 db Constant Bandwidth Curves for the Loaded and Null 
Adjusted Symmetrical Parallel-Tee Network 
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Fig. 13 Constant Q Loci and Locus of Symmetrical Frequency Response 
for the Loaded and Null Adjusted Symmetrical Parallel-Tee 
Network 
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11. RC Lattice Network - Complex Zero Synthesis and Sensitivity Analysis 
by Parameter Plane Techniques 
a. Introduction 
The RC lattice network which is investigated is shown in 
Figure 14. The parameters & and ‘ are real and variable in the arbitrarily 
selected range 0< & , 8 < 2.5. The transfer function of this network 


is given by: 


E (Pp) XB p” (oC eee peel (11-1) 


E; (p) x8 r+ (A&B +2+1)p+1 





where: p = RCs 
Frequency normalization makes the results which are obtained universally 
applicable. 

Parameter planes which are obtained graphically display the poles and 
zeros of (11-1) which occur for any choice of parameters & and 8 within 
the —— specified. These parameter planes show in what manner the poles 
and zeros change as the parameters are varied. Using the parameter planes 
one may select parameter values which are required in order to place poles 
Or zeros at desired p-plane locations within the capabilities of the net- 
work. The analysis performed shows that it is possible to readily place 
zeros anywhere in the complex region of the p-plane. 

Macroscopic root sensitivity can be determined by inspection of the 
parameter planes. Equations are obtained as described in Section 8 which 
permit microscopic root sensitivity to be computed for selected root loca- 


tions. 
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Be Poles and Zeros of the Transfer Function 
The poles of the transfer function are given by the roots 
of: 
ear + (AR +/8+1)p+140 (11-2) 
Since the network,ef Figure 14 digepameine RC, athe uoets of Gyal-2) all lie 
on the negative real axis in the p-plane. Let p= @ be such a root of 
(11-2). Making this substitution into (11-2) and rearranging the result- 


ing equation into the form of a conic section in the & -@ plane one obtains: 


(oF * +0) x - op + (O@ + 1) = 0 (11-3) 
Equation (11-3) represents a family of hyperbolas in the o& -B planes 
one hyperbola results for each value of @ . This equation is solved for 
A as a function of o& and @ and used to graph the curves displayed in 


Bagure 15: : 


@ - Cato ae (ile) 
ee i aes 


The zeros of the transfer £unctron are ene Loces OL. 


“8 p+ (a8 +1-f)p+1=0 (11-5) 
Following the procedures of Section 5 the constant zeta and omega loci are 
obtained in non-parametric form: 


(1) Constant Zeta Loci 





e(W) = 2 F XBW + (B- 978 -1+020 (11-6) 
g(W) = OW + (@-MB -1)W+2F =0 (11-7) 
From which: 
| 22 (p-.S-1) O O 
©) 235 48 (2 —vsel 0 
O . (11-8) 


Oo 
oo 
Lv 
t 
“ 
ou 
t 
}— 
ed 
tO 
oJ) 
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Expanding (11-8) one obtains: 
2(x , 6 5) =-4 5 8 (x 2.2041) 8 * + (-4 ba “eee 2) ort 1) reg aes) 


Thus: 
a= “ot = 2042) a CR Fl aed ay? - heal) (11-10) 


2( 7-20 +1) 
gives the constant zeta loci everywhere except where 43574 = Q. 
Equation (11-10) was used to obtain the constant zeta loci shown in Figure 
ie, 
(2) Constant Omega Loci 


Using Theorem II, Appendix II, for the cases j = 0 and j = 1, one 


obtains: 
-(8 - X68 -1)W + of 8 uw) Uy 
=0 Gist» 
(8 w*-1 ) 0 U, 
From which: 
WH 4B 4) ) = Kw (xAw 1) = 0 (11-12) 


Thus Xf 2 = 1/ a)? gives the constant omega loci everywhere except where 
of 8 SG = 0. Equation (11-12) was used to obtain the constant omega loci 
shown in Figure 16. 
oe Root Sensitivity 
(1) Macroscopic Sensitivity 
Root variations due to large parameter variations are read- 
ily observed by inspecting Figures 15 and 16. In the case of the sensitivity 
of the zeros of the transfer function, one can see that zeros with 0.3 <¥ 
< 0.7, and 0.5< W < 1.0 are relatively insensitive to parameter varia- 
tions. Given a specification on permissible zeta and omega variation, one 


can determine, practically by inspection, the required & and B8 tolerances 
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which will meet the specification. 
(2) Microscopic Sensitivity 
Zero sensitivity due to infinitesmal variations in the para- 
meters may be readily computed from formulas derived from the non-parametric 
equations of the constant zeta and omega loci as discussed in Section 8. 
Substituting derivatives obtained from (11-9) and (11-12) into (8-25) and 


(8-26) one obtains: 


Con 02 


Wd) (11-13) 
ya Diet 
OWi=- W (11-14) 
0B 286 


where: K(4aGW* + 2) #0 


OF 1 _ -F (axe -367 674406 7-448 - B42 B-1) (11-15) 
dk 2K (8B *wR - x? B72 +2~87% B*-2x8 428 -1) 


OF ; Fg gta -307 67 +648 7-3 B*408 446 -1) (11-16) 
af 2B XS -x? 67 4208 7- 87-2 WB 428 -1) 





where: f3 (8x & 8 -«7 8 242 o /S 12 oo wW8+28 -1) #0 
d. Interpretation of the Curves 
Inspection of the parameter plane for the zeros of the 
transfer function shows that it is possible to place a pair of complex con- 
jugate zeros anywhere in the s-plane by appropriately choosing the para- 


meter values & , & and the frequency normalizing factor RC. 





fe) CG @ 
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Fig. 14 RC Lattice Network with Variable Parameters 
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Fig. 16 Poles of the Variable Parameter Lattice Network 
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12. Application of the Method to other Problems 
This section contains examples of several natok s mhawe transfer 
functions are examined to determine the applicability of parameter plane 
methods. 
a. Pi Network 
The pi network shown in Figure 17 is useful for matching a 


resistive source to a resistive load at a fixed frequency. The source and 


load resistances are chosen as variable. 





Fig. 17 Pi Network with Variable Source and Load Resistances 


The transfer function of this network is given by: 





E(s) G 
a x3 reer + (XC UL +8C,L)s* + (x8 c, + o/3 Ci+ L)s +(% +8) 


(12-1) 
Equation (12-1) shows that the parameters appear in the coefficients in a 
manner that permits a parameter plane to be constructed which will show 
how the network poles are determined by the values of the variable para- 
meters. 
b. Unsymmetrical Parallel-T Network with Variable Parameters 
The network shown in Figure 18 is of particular interest here 
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because of the way that the parameters appear in the coefficients of the 


transfer function. 


RB wal 
@ 8 
@ C/X 
E; C/p AR E, 
Y EO 
p = RCs 


Fig. 18 Unsymmetrical Parallel-T Network with Variable Parameters 


The transfer function of this network is given by: 





Pol) = _ of @ pri t( eB? +6 "yp X18 + Bp + x8 (12-2) 
_~ x8 p> +( x8 +B “4X B+X+ B)p“+( A 3 2. Is a XB+%X+G)p t+ 1 


First inspection of (12-2) indicates that the pa’ terms will complicate the 
application of parameter plane techniques. The following is noted however: 
(1) The network is passive RC, therefore the poles are all real and lie on 
the negative real axis in the p-plane. Hence p = @~ may be substituted into 
the denominator of (12-2). This denominator is then equated to zero and 

the resulting equation solved by the quadratic formula for R as a function 
of eo and & . Thus the parameter plane for the poles may be constructed. 
(2) The quantity ®8 may be factored from the numerator of (12-2). The 
resulting polynomial is of the proper form to be analyzed by parameter plane 


methods. 


(i, 





Cu Bridged-T Network with Variable Parameters 


The bridged-t network shown in Figure 19 can be used to gener- 


ate complex conjugate zeros in the left half s-plane. 


BR 





Fig. 19 Bridged-T Network with Variable Parameters 


The transfer function of this network is: 


ee) = was) ar + (M+ 1)s + 1 (12-3) 


E, () od) Ea + (X+6+ Dig 2 





where p = RCs 
The form of Equation (12-3) readily permits parameter plane mien and 
design. 
d. Bridged Twin-T Network with Variable Parameters 
The network shown in Figure 20 is an RC Bridged Twin-T 


Network with variable bridging and shunting elements. 


x A i 
e = ‘A O 


IE 3 


SR YC 


Fig. 20 Bridged Twin-T Network with Veriable Parameters 
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The transfer function of this network is given by (12-4): 


ae) “> (12-4) 
E, (p) 


Br) Pt (26BS+ 26YS 5) PHYO 20aSHBx'+ 25 )p%4(2e10/ar25) Pt! 
(Ker'S) PT (2d ar5 +2BVS +2485 +4Ay' +¥'S) P?+ 
(24BF +A y¥'+208 t4Bdr2 YS +8 EY 42d) poe 
(X38 +28+25+2)P4/ 





Where p = RCs 

Inspection of (12-4) reveals that one may construct several parameter planes 
depending on which pair of the quantities o& , S > yy , db. are considered 
variable and which pair are considered fixed. The following possibilities 


are available: 


Variable Parameters Fixed Parameters 
Xo , 8 ard, 
ie a, 
X,Y" p,o 
prt o, ¥ 
Lod et 
B>¥ Ad 


e. Capacitive Divider Matching Network 


The network shown in Figure 21 is frequently used for 
frequency selective coupling between high output resistance and low input 


resistance amplifier stages. 





Fig. 21 Capacitive Divider Matching Network 
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The transfer function of this network is: 


Ne - Bx wee (12-5) 
Es) (BYd Lye BOS 17 +L) s°+( fr LL iee+y@e 1 RY Ore 





(xy! L+ ABD + WY + L)s +o 
Equation (12-5) shows that it is possible to construct parameter planes 
with the capacitances ¥ and @ variable and the source and load resistances 
A and = fixed or vice versa. 
iy Two Section Resistively Loaded LC Low Pass Filter 
The LC low pass filter shown in Figure 22 has variable 


source and load resistances. 


Fig. 22 Two Section Resistively Loaded LC Low Pass Filter 


The transfer function of this network is: 


EGS) 9. 6 (12-6) 








E_(s) 4 3 
Oo (B L,L,C,C,)s + (6 CCL +L LC ,)s + 
Zz 
of + + 
(BL C+ BLC)+ BL,C,)s° + (98 C+ Bc, L,tL,)s 


(+B ) 


The form of this function permits application of parameter plane techniques. 
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g. Butterworth Filter Pole Sensitivity to Coefficient Variations 


Parameter plane techniques can be used to investigate the 
effects upon pole locations of Butterworth Filters due to variations in the 
coefficients of the transfer function caused by component tolerances, aging, 
etc. As an example, consider the transfer function of a normalized 4th 


order Butterworth filter: 


Bo(P) : (12-7) 





Bir) SE Ned tains 2. Gisocd 
Assume that the last two coefficients in the denominator of (12-7) are 
variable. Denote these variable coefficients as o& and 8 and equate the 


denominator of (12-7) to zero. This results in the polynomial equation: 


at + 2.613p> + gees of pt+f = 0 (12-8) 
The parameter plane for this polynomial is displayed in Figure 23. When 
OX = 2.613 and fs = 1.0 as is the normal case, the poles of (12-7) are 


equally spaced on the unit semicircle in the left half p-plane. By enter- 


. 
oo 


ing Figure 23 with values of & and ik , the resulting poles of (12-7) 


may be read directly. 
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13. Algebraic Design of Electric Networks by Parameter Plane Techniques 
a. Introduction 

Parameter values which will guarantee that certain poles and zeros 
of a rational network function appear at desired s-plane locations and at 
the same time guarantee that other specifications are met may be algebraical- 
ly determined by employing parameter plane techniques. This may be accomplish- 
ed without the necessity of drawing parameter plane curves and thus, more 
than two adjustable parameters may be simultaneously considered. In general 
it is necessary to have at least n adjustable parameters in order to make n 
specifications. While the design technique does not mathematically demand 
it, some preliminary investigation is desirable to determine whether or not 
a physically realizable set of parameter values is likely to be found which 
will meet the specifications. If such an investigation is not performed 
then it may be the case that the design technique will generate only sets 
of parameter values which are unrealizable. 

eae Types of Specifications Which May be Used 

Any specification on the network function or on the network be- 
havior which can be expressed as an algebraic equation in the adjustable 
parameters can be used in conjunction with this technique of design. These 
specifications can include both pole-zero locations and performance specifi- 
cations. Some of these are: 

(a) Pole-zero locations. The parameter plane equations of Section 
4 can be used to specify desired pole-zero locations. 

(b) Zeta of a pair of complex roots may ‘be specified without 
specifying the corresponding omega, or vice versa, by using the non-para- 
metric equations of the constant zeta and omega loci presented in Section 


a 
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(c) Bandwidth. The equations of the constant bandwidth curves 
developed in Section 6 can be used to specify bandwidth. 

(d) Q. The equations of the constant Q loci, Section 7, may 
be used to specify Q for second order systems. 

(e) Maxima and minima in the frequency response can be specified 
by equating the derivative of the magnitude function to zero at the fre- 
quency where the maximum or minimum is desired. 

(f) Infinite frequency or DC gain of a network may be specified 
by expressing these quantities as algebraic equations in the adjustable 
parameters. 

Cc. The Technique of Algebraic Design 
(1) Pole Specifications 
Consider either the numerator or denominator of a rational 


network function written as the polynomial equation: 


3 as) = 0 Mis.) 
k=0 
where the coefficients a, are functions of the adjustable parameters q,: 


ae = £ Cd) 9doo+++945) (13-2) 


Parameter plane theory provides a number of algebraic equations which 
relate the root factors Pw, eo to the parameter values q;° These are: 
(a) The basic parameter plane equations for real and 
complex roots: 
1. For real roots: 


ao set. (13-3) 


x 
iM 
oO 


2. For complex roots: 
k k 
yi (-1)"a, w Up 64 ) =0 (13-4) 
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where j is any integer. 
(b) The constant zeta and omega loci in non-parametric form 
as developed in Section 5: 


z(a.,% ) = 0 (13-5) 


W(a,, w »}=0 (13-6) 


(2) Other Specifications 


Consider the rational network function (13-7) which may be 


the transfer function of some network: 


(13-7) 
k 


ry bs 
k=0 


where the coefficients ay and be are different functions of the adjustable 


parameters q,: 


a = £ Cd) 9dos-+29d5) (13-8) 
by = 8.04, 9dos-+- 4.) (13-9) 


Specifications on the behavior of F(s) must be written as algebraic equations 
in the parameters qs in order to apply this design technique. This is 
straightforward for the common types of specifications. Sections 6 and 7 
show how this is accomplished for specifications of Q and bandwidth. 

(3) Solution for the Parameter Values which Meet Specifications 

If r poles and zeros of the network function are specified, 

then Equations (13-3) through (13-6) may be employed as necessary to gener- 
ate r independent equations in the adjustable parameters qs: These r equa- 
tions taken together with perhaps s other equations arising from other 
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specifications form a set of (str) equations which must simultaneously be 
satisfied. It is presumed that the number of adjustable parameters q; 

is equal to or exceeds the number (str) of specifications. In order to 
obtain discrete solutions for the parameter values, the number of specifi- 
cations must be made equal to the number of adjustable parameters. This 
may be achieved by fixing as many of the adjustable parameters as required 
to obtain this equality. Once the number of adjustable parameters equals 
the number (str) of specifications, the set of (str) algebraic equations 
may be simultaneously solved. Simple cases may be solved by hand, however 
the equations are generally complicated and require computer solution. 
Several iterative procedures are available to achieve this, and are usual- 
ly some type of search routine which seeks parameter values which will force 
all equations to zero simultaneously. 

In general the set of (str) equations is non-linear and its solution 
will result in several sets of parameter values which mathematically satisfy 
all equations simultaneously. Some of these solutions will be rejected be- 
cause they are unrealizable. It should be kept in mind that the mathemati- 
cal solution may not always be physically acceptable even though realizable. 
This situation may occur, for example, when specifying a maximum in the 
frequency response characteristic by selecting parameter values which will 
cause the derivative of the magnitude function to be zero at some specified 
frequency. The mathematics may indicate a solution which results in a 
minimum rather than the desired maximum. In any event solutions should be 
checked to ensure that they meet all of the specifications. If a solution 
exists which is both realizable and physically meets the specifications, 


then it will be found among the solutions to the (s+r) equations. 
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d. Example 


Consider the case of the symmetrical parallel tee network 
which was investigated in Section 9. The zeros of the transfer function of 


this network are given by the roots of: 


o /3 p + QE eh et sO 


Since this equation has two variable parameters, two of its three roots may 
be specified. A preliminary analysis has indicated that complex roots are 
possible with real and positive parameter values. Assume that the para- 
meter values are desired which will cause a pair of complex roots to have 


zeta equal to 0.2 and omega equal to 1.0. The basic parameter plane equa- 


tions for complex roots are chosen with j O and j = 1. These become: 


0.84 ¢%8 -1.2 8 =0 


1.2 A = 0.84 


Solution of these equations produce the desired parameter values: 


=i 
8 


By dividing the quadratic factor associated with the specified complex roots 


1.44 


0.70 


into the cubic equation (13-10), the real root is determined to be located 
at p = -l. 

ec. Example 

Consider again the parallel-tee network of the previous example. 

Instead of completely specifying the complex zeros, let only the omega be 
specified, and instead of specifying zeta let the value of the real root be 
specified. Assume that the omega of the complex zeros is specified as 1.2 
and that the real zero is to occur at s = -0.5. In order to obtain the 
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desired parameter values which will result in these roots, the equations of 


the real roots and the constant omega loci are employed: 


= ao F = 0 


k=0 
Wa, WwW) = 0 
From Section 9, and with @ = 1.2 and T = -0.5, these equations are: 
-0.125%8 + 0.508 - 1.08 + 1.0 = 0 (13-11) 
2 2 2 
2.99 4° 4@~ - 4.14 x8 + 2.88 8 - 1.0 = 0 (13-12) 


Equation (13-11) is solved for & as a function of f° : 
a = 0.508 - 1.0 (13-13) 
=O2l2a 73 
Equation (13-13) is substituted into Equation (13-12) to form the quadratic 
in B: 
64.3387 -221.34 + 190.10 = 0 (13-14) 
Equation (13-14) is solved for two values of & and then Equation (13-13) 


is used to obtain the corresponding values of % . The results are: 


xe £ : 
1%) 0.85 GB, 1.65 


oY, = 0.46 (3, = 1.79 


Upon checking of these results a, and B, are found to produce the desired 
zeros. A, and 3 5 are rejected since they produce three real zeros, one 

of which has the desired sigma, rather than one real root and a pair of com- 
plex conjugate roots as desired. This situation was anticipated as discuss- 


ed in the introduction to this section. 
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ae Example 


In the previous two examples only zeros of the transfer func- 
tion were specified. In this example a pole is specified and a constraint 
is placed on the zeros. The symmetrical parallel tee network discussed in 


Section 9 is used as a vehicle. The transfer function of this network is: 


E(P) _ 6p +2B8p+2@p+r (13-15) 


E (p) Bp + (2R+2%/B4+X )pa(2 B +e + 2)ptl 
where p = RCs 
Assume that a pole is required to be placed at p = -0.30 and that a pair of 
complex conjugate zeros is to have a damping ratio rf = 0.4999. The equa- 
tion of the constant zeta loci for the numerator polynomial which was obtain- 


ed in Section 9 is used to express the constraint on zeta: 


((8.00X -16.00) & +8.00)8 = 0 (13-16) 


The equation for the constant sigma curves for the denominator polynomial, 


also obtained in Section 9, is used to express the constraint on the pole: 
- 0.153% -0.420 -0.050&% +0.400 = 0 (13-17) 


Equation (13-16) is solved for & as a function of @ : 


of = 2.00 - (1.00/86 ) (13-18) 


When Equation (13-18) is substituted into Equation (13-17) the following 


quadratic results: 


Q* + 1.528 - 1.84 = 0 (13-19) 


Equations (13-18) and (13-19) give the following pair of parameter values 
which mathematically satisfy the constraints placed on the pole and zeros 


of the transfer function: 


0.74 3B. = 0.80 


a 1 


Xx -2.31 


| 
No 
rr 
Ww 
ms 
No 
ll 


Z 





Upon checking, oT and 8 , are found to satisfy the specified constraints. 
x, and (2, are rejected as unrealizable since (O 9 is negative. 
Qe Example 

Consider the loaded and null-adjusted symmetrical parallel 
tee network discussed. Presume that one wishes to select a value of the 
source resistance & , and the load resistance s , which will cause a 
pole to be located at p= Q = -0.3 and which will result in a value of 
Q =0.34. Referring to Section 10, Equation (10-2), it is seen that the 


specification on the real pole requires that: 


-0.84 a B - 0.11 & -0.11 8 + 1.40 = 0 (13-20) 
Equation (10-10) of the same section gives the equation which must be satis- 


fied to meet the specification on Q: 
aus 2 2 2 2 
1.8496 X" B™ - 0.3008 XA~ + 0.8496 BX - 0.3008 %' AB - 4.4512 xB - 


0.15048 + 0.8496%~ - 0.1504~ + 0.4624 = 0 (13-21) 
Equations (13-20) and (13-21), when simultaneously solved by digital computer 
produce the solution: 
oe = 1.95 
B = 0.677 
These values are physically realizable and meet the specifications. 
h. Conclusions 
The parameter plane technique presented in this section is 
the only design technique known to the author which allows parameter values 
to be determined which will cause both pole-zero and performance specifica- 
tions to be simultaneously satisfied. The types of specifications which 
may be made are limited only by the requirement that they be expressible as 
algebraic equations in the adjustable parameters. For common specifica- 


tions this is readily achieved. 
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14. Conclusions and Recommendations for Further Work 
a. Form of the Coefficients 
The work which has been presented in Section 4 shows how parameter 
plane equations have been developed for polynomials with coefficients of 


the form: 
a, = be + C. a + d x6 + fo (14-1) 
This extension of the theory permits a wide class of network functions to 
be analyzed. There remains, however, cases where higher order terms in 
and appear in the coefficients. These include terms Mieeme AP alee ; 
a 3 g etc. Work currently being performed by Lt. A. R. Miller, USN at 


the U. S. Naval Postgraduate School has shown that it is possible to solve 


the parameter plane equations for coefficients of the form: 
=b,X~+e 97 +4 %8 +e +f B+ (14-2) 
*;, ak ke Ok k k ae 
Preliminary work done by Lt. R. E. Hudson, USN, also at the U. S. Naval 


Postgraduate School, has shown that it is theoretically possible to solve 


the parameter plane equations for coefficients of the form: 


a= bX" 6 7 + C48 +d. x /3 oa ex. + £ B° + go + hot qa, 6 + t 
(14-3) 
There remains however much to be done, particularily in interpreting the 
curves which result from solution of the parameter plane @quations associated 
with (14-2) and (14-3). 
It is, of course, desirable to allow even higher order combinations of 
o and & to appear in the polynomial coefficients than allowed by (14-3). 


At the current state of development it is not clear just how to achieve this. 
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b. Non-parametric Equations of the Constant Zeta and Omega Loci 


The derivation of the non-parametric equations of the constant 
zeta and omega loci presented in Section 5 permits these loci to be more 
readily constructed than was previously possible. More important, however, 
it forms the bridge between parameter plane theory and the theory of plane 
algebraic curves. To date this relationship has not been seriously investi- 
gated. It is the author's opinion that such investigation will lead to new 
and valuable information about the behavior of the loci. Such information 
may well lead to the development of sketching techniques which would con- 
siderably simplify the construction of the loci. Additionally it might 
permit new insight into the behavior of the zeta equals zero locus which is 
the stability limit for the poles of active networks. 

Cc. Constant Bandwidth Curves 

Simultaneous correlation between frequency response character- 

istics and root locations upon the parameter plane is now possible as a 
result of the development of the constant bandwidth curves which were develop- 
ed in Section 6. Families of such curves permit the frequency response of a 
network to be sketched for any choice of parameter values within the range 
for which the curves are drawn. 

The equations developed in Section 6 apply to coefficients in the trans- 
fer function of the form (14-1). When parameter plane techniques are 
developed for coefficients of the form (14-3) it will be necessary to extend 
the work of Section 6 to include this coefficient form. Preliminary work 
done by the author has shown that this is possible by following the same 


procedure as employed in Section 6. 
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di. Constant Loci 
As defined by Morris (4), Q is a mathematical parameter of a 
second order characteristic equation. Thus defined, Q may be employed to 
characterize RC networks for which the conventional definitions fail to 
provide useable correlation between the network characteristics and the 
value of Q. 

Loci of constant Q, developed in Section 7, permit the value of Q 
associated with any choice of network parameters to be read directly from 
the parameter plane. Alternately one may select parameter values from the 
parameter plane which will achieve a desired value of Q within the capabili- 

ties of the network. 

é. Algebraic Design by Parameter Plane Methods 

The algebraic design techniques presented in Section 13 per- 

mit m network parameter values to be determined which will guarantee that m 
specifications upon the network are met. These specifications can be either 
pole-zero locations or other specifications which can be written as algebraic 
equations in the m adjustable parameters, or both. As pointed out in 
Section 13, it is necessary to have the same number of adjustable parameters 
as specifications in order to obtain discrete sets of parameter values as 
solutions to the problem. Unless these specifications are realistic ones 
in terms of the network capabilities, then it may be the case that no 
physically realizable solutions will be produced by the design technique. 
It is more likely that a realizable solution would be found if the system 
had more adjustable parameters than specifications and if the specifications 
were stated as inequalities rather than strict equalities. For example, one 
might specify Q greater than 0.34 rather than strictly equal to 0.34 provid- 
ing that this is acceptable. This procedure would lead to a set of equations 
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which would contain some strict equalities and other inequalities and hence 
other methods than those of Section 13 would have to be applied in order to 
obtain parameter values which would meet the specifications. The author 
suggests that non-linear programming techniques be investigated in connec- 
tion with the solution of such a set of equations. 

1A Summary 

The basic parameter plane techniques, and the extensions to 

them presented here, provide a powerful tool for the analysis and design of 
variable parameter, active and passive, electric networks. These techniques 
are capable of providing information about the poles and zeros and sensitivity 
of variable parameter networks which is not possible by any other method 
currently known. It is, for two parameter problems, all and more that 
Evan's root locus technique is for single parameter problems. Like the 
root locus technique it is applicable to a much wider sphere of problems 
than that of electric networks. Subject to the restrictions on the form of 
the coefficients previously discussed, parameter plane techniques may be 
used to investigate the singularities of any variable coefficient rational 


function in the complex variable s. 
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APPENDIX I 


Inductive proof that: 


UC FU, CF ) = CF 8 7 ee D (1-1) 
for all integers i and }j. 
The inductive proof that (I-1) is true consists of the following steps: 

1. Show that (I-1) is true for i = 2 and i = 3. 


2. Assume that (I-1) is true for i = m-l and i =m. 


3. Show that (I-1) is true for i=m+H1, given that it is true 

for i = m-l and i = m-2. 
Step la: 

1 Siz 

U,(FU, CF) - U,CF ULF ) = -U,_ CF ) (1-2) 
but 

U(%) = 27 

U,¢ f mS 1 

Therefore (I-2) may be written: 

27 Ure Fg) 0,0 FS WCF D (1-3) 


but UL Z ) = -U, ¢ Zz ) 

Hence (I-3) becomes: 

2 FB ing) =, CFO) 58 e F (1-4) 

(I-4) is recognized as the recursion relation which defines the 
Chebyshev functions. Thus (I-2) is true for all j. 

Step lb: 

i= 3 
U3(4) UC) - UWF WF ) = -Us_ C5) (1-5) 
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but : 


UF) = 4 5 7-1 


Therefore (1-5) may be written: 
2 
GS-1, CF) - 27 uF) =u, CF) (1-6) 


Making use of the recursion relation for the Chebyshev functions, (I-4), 
it is possible to write (I-6); 

(AFP, CF - 4 FU CFI + 2s CF) =U; 3%) 1-7) 
Expanding (I-7) and cancelling where possible), produces the recursion relation 
for the Chebyshev functions, (I-4). 

Therefore (1-5) is true. 

Step 2: 


Assume that the following equations are ture: 


Ba 5 Us OF) - Gyo SD BD ee ea F? (1-8) 
UF) UF) - UL CF WCF) = UF) (1-9) 
Step 3: 

Un- p41 5 Un- gti ¥ ) (I-10) 


Applying the recursion relation for the Chebyshev functions, (1-4), to the 


left hand side of (I-10) and rewriting, one obtains: 
2 Flu, CF) - CH) =U sd (I-11) 


Substituting (I-8) and (I-9) into (I-11) one obtains: 
25 FU Fae ee UC ZI - (U1 (FU, _4(F)-U,_9( FU, (3 J = 


Une H1( 3) | (I-12) 


Equation (I-12) may be written: 
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(250 (FO (FM, CF) “125 Oe CF -W CF)T OCF) = BF) 
(1-13) 


By use of the recursion relation for the Chebyshev functions, (1-4), Equa- 
tion (I-13) may be written: 
Be O78 > Bal GOGO) SS) 


Thus (I-1) is true for i = mtl, given that it is true for i = m-l and i = 


(I-14) 


m-2. This concludes the proof that (I-1) is true for all i and j integer. 
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APPENDIX II 


Proof that: 
n k k 
a (-1) a, w UF) = 0 (II-1) 
k=0 
for all integers j and k 
Theorem: 
= k ok 
If: 2 (-1)"a,wu( 4s) = 0 (II-2) 
k=0 
and if: 
n k k 
> (-1) "aw U, 403 ) =0 (II-3) 
k=0 
Then: 
RY '(.1)*a UC) Om tor eel) eaeoeee 
5. . k= j i) all integers 
saad j and k (11-4) 
Proof: 


Subtract (II-2) multiplied by U;_4¢ Z ) from (II-3) multiplied by 


UC 3 ). This gives: 


"\ 
>, CD*a,w* ULF U, (3) = Y Fw CF] = 0 (11-5) 
k=0 


Employing the results obtained in Appendix I, namely (1-1), Equation (II-5) 


may be written: 


Nn 
>, Cv*aak (0, (5) = 0 (II-6) 
k=0 


Multiplying (II-6) by -1 concludes the proof of (II-4). 
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APPENDIX III 
Table of the Chebyshev Functions 
The Chebyshev Functions TC "i ) and U,( z ) are defined by the following 


recursion relations: 


Tea ) -2 FHCF) + HOF) #0 


Te? a! 
17 aon) 
Ver F > -BUCF +07) = 0 
Uae 0 
hy) al 


The first few of these functions are tabulated below: 


k Tas) U.(F ) 
0 l 0 
1 i] 
2 2a at Ars 

3 Z 
3 45° -3% 43 -1] 
4 8 y*-8 541 87° -4% 
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